Dynamical formation of floating Rydberg excitation crystals 
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The dynamics of a cloud of ultra-cold Rydberg atoms is studied at off-resonant laser driving. We 
find that excitation crystals are formed dynamically as a consequence of interaction-induced resonant 
excitations. These crystals can extend over the entire ensemble volume, have lattice constants 
independent of the ensemble length, are not localized relative to the ensemble, and sensitively 
depend on the shape of the interaction potential. Compared to previously proposed crystals, these 
distinct properties lead to qualitatively different predictions for the spatial excitation density, the 
Mandel Q parameter, and the total number of excitations. 
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During the last few years an immense research activity 
both experimentally and theoretically aimed at a detailed 
understanding of the properties of ultra-cold gases whose 
atoms are excited to states of high principal quantum 
number n. Such highly excited Rydberg atoms have ex- 
treme properties which make them promisiiig candidates 
for a number of fascinating applications [l|, |2[ . Most im- 
portantly, they feature long-range van der Waals (vdW) 
interactions, which lead to interesting effects like the 
dipole blockade [3|, |4|. After initial work mainly con- 
cerned with bulk properties of Rydberg ensembles, cur- 
rent effort focuses more and more on spatially resolved 
observations. This is motivated not least by experimen- 
tal progress [5] and theoretical proposals [a, lZ| for spa- 
tially resolved excitation imaging. Also on the theoreti- 
cal side, in particular rate equation models [8[ |9| and full 
many-body simulations on truncated Hilbert spaces [lOl - 
14| provide a handle to access spatially resolved proper- 
ties. 

In this spirit, it has been predicted that spatial pair 
correlations can be induced in a three-level Rydberg gas 
via the so-called anti-blockade arising from resonant ex- 
citations due to single-atom Autlcr-Townes splitting [15J . 
These correlations could be measured based on mechan- 
ical forces due to the vdW interaction. The induced 
particle motion leads to an encoding of position cor- 
relations into a time-dependent Penning ionization sig- 
nal 16]. Interestingly, in this way, spatial information is 
gained without a spatially resolved measurement. Fol- 
lowing these works on spatial properties of pairs of Ry- 
dberg atoms in large ensembles, it was subsequently 
shown that also the Rydberg ensemble as a whole could 
form cry stalline structures in the Rydberg excitation den- 



sity [17|, |18|. In a repulsive vdW gas of Rydberg atoms. 



it can be energetically favorable for a given number of 
excitations to assume a highly ordered, crystalline state, 
with distances between the excitations maximized. If the 



laser detuning acts opposite to the interaction contribu- 
tion, these ordered states become the quantum mechani- 
cal ground state of the system. It was proposed that such 
ground-state crystals (GSC) could be produced using 
chirped laser pulses [ITI . Il8| . The chirp of the laser driv- 
ing induces adiabatic passage to the energetically most 
favorable state corresponding to the crystalline excita- 
tion structure. Recently, it has been shown that spatial 
correlations also arise for resonantly driven atoms, which 
can be understood from entropic arguments |19j . 

Here, we show that Rydberg excitation crystals ex- 
tending over the entire ensemble volume can also arise 
dynamically for ofF-resonant laser excitation, without the 
need for laser chirping. For this, we study in detail the 
spatially resolved properties of a one-dimensional disor- 
dered gas of Rydberg atoms via numerical simulations 
on a truncated Hilbert space. The crystalline structures 
predicted here are fundamentally different from the GSC 
reported previously. First, their lattice constant does 
not depend on the trap geometry. Second, the crystals 
are not spatially localized with respect to the trap, but 
can float over a certain position range between differ- 
ent realizations. Third, common observables yield re- 
sults qualitatively different from those for the GSC. In 
particular, due to the floating nature of the crystal, the 
total and the excitation-number resolved spatial excita- 
tion densities exhibit characteristic structures. Methods 
to experimentally verify the crystalline structure are dis- 
cussed. The underlying resonant excitation mechanism 
is dynamical, inviting time-resolved studies of the forma- 
tion of spatial correlations. Part of our predictions can 
already be probed through the detuning and density de- 
pendence of the bulk total number of excitations without 
the need for spatial resolution. 

Our model system is a one-dimensional cloud of 
Rydberg atoms in two-level and frozen-gas approxima- 
tions (2[. We assume temporally and spatially constant 



laser intensity and wavelength. The corresponding many- 
body Hamiltonian in a suitable interaction picture and 
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rotating wave approximation reads {h= 1) [3, llCllll[ 
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with s]^'p = |a)i (/3| and cri*^ = sig 
first part of this Hamiltonian contains a sum over single- 
atom contributions. It includes the detuning A between 
laser frequency and atomic transition frequency, and the 
laser coupling between the ground state \g) and Ryd- 
berg state |e) with Rabi frequency 51. The second part 
accounts for the vdW interactions between two atoms 
in the Rydberg state at a mutual distance rij. We nu- 
merically solve the time dependent Schrodinger equation 
with Hamiltonian ([Ij for given atom positions r^ and pa- 
rameters Cg and A, starting from an initial state with 
all atoms in the ground state. For a many-body sys- 
tem, a naive integration is impractical, as the number of 
states grows exponentially with the particle number. We 
overcome this problem by exploiting the Rydberg block- 
ade to truncate the Hilbert space to a physically rele- 
vant subspace 11]. All simulation results are checked 
for convergence with respect to the truncation parame- 
ters. Additionally, to simulate larger ensembles and to 
obtain a convergence of the ensemble state with time, 
a Monte Carlo sampling over runs with different atom 
positions is employed [ll|, |20|, 



2l|. In all the simulations 



shown here, the Rabi frequency was kept at f2 = 10 MHz, 
and the interaction strength at Cg = 6.6THz/im^. This 
choice, however is not restrictive since the results can be 
mapped to any other parameter values by rescaling time 
and length. A^ = 30 atoms over a length L = 30 /Ltm corre- 
spond to a line density oi p = 10'' cm~^, which translates 
into a volume density of lO^^cm"^. For these param- 
eters, the systems studied in this work would typically 
converge at physical evolution times of about 2 fis, while 
the reported observables have been obtained at time 5 /xs. 
Our main observable is the pair correlation function 
g^'^' . It is a measure for the conditioned probability for 
an excitation if another excitation is already present at 
a distance r. It can be calculated as |22| 
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where ^^ sums over pairs with distance in [r, r + Ar] 
and Np{r) is the total number of such pairs. 

We now turn to the results. Starting point of our 
analysis is the observation that pronounced resonances 
emerge in g^^' at regularly spaced distances for detunings 
A > J7, see Fig. [TJ In the following, we show that these 
resonances can be traced back to a crystalline excitation 
structure. At large detunings (A/J7 ^ 1) and without 
the interaction (Cg — >■ 0), the laser field is off-resonant 
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FIG. 1. (Color online) (a) Pair correlation function g^ '{r) as 
a function of inter-particle distance r and detuning A. Dashed 
lines are theoretical predictions for resonances explained in 
the text, (b) Sections through (a) at detunings A/Q = 3, 5, 7. 
Parameters are L = 30 fim, N = 30, fi = 10 MHz, and Ce = 
G.eTRzfim^. 



with all transitions and does not induce excitations. But 
if the vdW interaction is included, it can shift transi- 
tions to higher excited states into resonance with the 
laser 15|, Il6l . l23j . From such resonance conditions, dis- 
tances r^es can be deduced at which atoms are resonantly 
excited for a given detuning. For example, the ground 
state is resonantly coupled to a doubly excited state via 



two-photon excitation if 2A = C'g/rJ?;,^^- This condition 
is shown as the lowest dashed line in Fig. [ija) coinciding 
exactly with the first maximum oi g^^' from the numerical 
simulation results for large A. We have identified the sec- 
ond maximum indicated by the next dashed line towards 
slightly higher r as the resonant transition from an m-fold 
to an m-t-1-fold excited state, satisfying A = Ce/rf^^i. 
If a pair is further excited to a triply excited state, the 
distance between atoms 1 and 3 is rres,2 + ?'rcs,i, causing 
the third resonance line in Fig. [TJa). Resonances due to 
higher excited states can be explained analogously. We 
can thus directly trace back the emergence of spatial cor- 
relations to resonant excitation channels. 

Next, we analyze the structure of the spatial corre- 
lations at high detunings in more detail. Figure [IJb) 
shows that the pair correlation resonances become nar- 
rower with increasing detuning. The origin of this effect 
is illustrated in the inset of Fig. [H^a). The resonances 
are of Lorentzian shape in energy space as a function 
of the detuning. The interaction potential V{r) trans- 
lates this dependence on the detuning into a distance 
dependence. For the same energy-resonance width, the 
position-space resonances become narrower as the dis- 
tance decreases at which the laser moves into resonance. 
Hence, an increase of the detuning leads to a reduction 
of the resonant width. For the particular case of a two- 
photon resonance, the energetic resonance width depends 
on the detuning itself. To quantify this, we calculate 
the steady-state value of the probability for double exci- 
tation Pee from the two-atom master equation. In the 
limit of vanishing spontaneous decay rate, this yields 




") 3 


/iT'"--.^ 


2.5 


/ \ X 


■a 

Z" 2 


X'' 


1.5 


/ ''^ 


1 


v_^.'''' 



1 

0.5 



. 



Ot 



-0.5 



2 4 6 
AJQ. 



FIG. 2. (Color online) (a) Dependence of the two-photon ex- 
citation linewidth on the detuning. The solid line is obtained 
from Eq. ([3]), red crosses are obtained from the numerical 
simulations. The inset shows the detuning-dependent trans- 
formation of energetic resonances into spatial resonances via 
the interaction potential V{r) = Ce/r^- (b) Total number 
of excited atoms and Mandel Q parameter as a function of 
the detuning. Solid lines show A^ = 30; dashed lines show 
A'' — 45. Cc and L are as in Fig. [U 
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FIG. 3. (Color online) (a) Rydberg excitation density in ex- 
citations per ura vs. position on the string and detuning. 
Dashed line: resonance line corresponding to 2 A — Ce/r^ 
as explained in the text, (b) Excitation density for A/fi = 7 
split into excitation- number subspaces. 



Pee = f^'*/{4[f74 + {V - 2A)2(f72 + A^)]}. Linearizing 
the potential around ri-es,2 we obtain a Lorentz curve in 
5r = r — Tics, 2 with the position-space width 



r = 



' res, 2 



3CeVn^ + A2 
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V02 + A2 
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As shown in Fig.[2lja), the numerical simulation data fully 
agrees with Eq. ^ up to an overall prefactor of about 
2 due to additional atoms in the trap which broaden the 
resonance compared to the idealized two-atom case. 

In summary, we find that in the limit of large A, 
the position-space resonance width decreases as A~^^/^. 
Therefore, in our model, resonant excitations are only 
possible at definite positions with lattice spacings deter- 
mined by the shape of the interaction potential. This 
leads to crystalline structures independent of the trap 
length, in contrast to GSC which have lattice constants 
determined by the trap geometry. 



After having established the formation of stronger spa- 
tial ordering with increasing detuning, we now address 
the question how this ordering manifests itself in vari- 
ous observables. The distribution of the number of ex- 
citations is characterized by the Mandel Q parameter 
Q = ((/V^d) - (A^yd)')/(^yd) - 1 [221. For a Poisso- 
nian distribution of excitation numbers, like in a coherent 
state, the Q parameter is zero, while for sub-Poissonian 
statistics it is negative. In previous proposals on GSC, 
for a given detuning, a definite number of excitations 
was predicted in the system. Then, ideally Q reaches the 
value —1 for a Pock state. In contrast. Fig. [2l^b) shows 
that in our setup, the Q parameter assumes a minimum 
at some positive detuning, but then increases again to- 
wards larger A. This is due to the fact that with increas- 
ing detuning, the spatial excitation resonances become 
narrower, as shown in Fig. [2ja). Consequently, the to- 
tal number of excited atoms decreases with increasing A 
since the distance range over which resonant pair excita- 
tion is possible shrinks. For higher densities, this effect 
is expected to set in at higher detunings, since at a given 
resonance width, it is then more probable for an atom to 
have a resonant partner at a suitable distance. This is 
supported by our numerical simulation data, as the maxi- 
mum of the number of Rydberg excitations A^ryd is shifted 
to higher detunings for higher densities, see Fig.IlJb). We 
have found that at large detunings, the system evolves 
into a bimodal excitation distribution, with one fraction 
in the ground state, and the second fraction distributed 
around a parameter-dependent excitation number. As a 
direct consequence of this bimodal distribution, the Q 
parameter deviates from the Pock-state result, and thus 
is not a definite indicator for the presence of spatial order 
in the case of resonant excitations. 

It is rather the spatial order itself which can be used to 
verify the crystalline structures. In the crystalline state, 
the distance between excitations predominantly is a mul- 
tiple of the resonant distance. But interestingly, this does 
not translate into a peaked structure in the spatially re- 
solved excitation density. Fig.[3l^a) shows that at high de- 
tuning, the excitation density has a rather complicated 
step-like dependence on the position. This can be ex- 
plained with a geometric argument: If we consider, e.g., 
the doubly excited states, we notice that these are only 
populated if the distance between the excited atoms is 
the resonant one. If we now ask how such a pair of exci- 
tations with fixed distance can be placed on the string, 
and assume that each of these possibilities is realized with 
equal probability, it becomes clear that atoms located 
within one resonant radius from the ends of the trap are 
excited only half as often as atoms at the center. This ex- 
plains the outermost edges in the excitation density. The 
dashed line in Fig. |3l^a) indicates the position ri.es,2 away 
from one edge of the string coinciding with the position 
of the step in excitation density, clearly supporting this 
interpretation. The other edges can be explained analo- 
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FIG. 4. (Color online) Illustration of the crystal formation: 
Positions of the excitations are determined from the exact 
quantum state by a Monte Carlo procedure mimicking a mea- 
surement. The distances of the excitations to the leftmost 
excitation are collected in a histogram. Note that the scale 
of the ordinate in (c) and (d) differs from (a) and (b). The 
height of the two-photon resonance is 12 in (c) and 19 in (d), 
respectively. Parameters are as in Fig. [T] A/Jl — 6. In (d) 
the steady state has been reached. Excitation times are given 
in fis. 



gously, taking into account higher excited states. This is 
also illustrated in Fig. ^h) in which the excitation den- 
sity is split up into contributions of states with different 
excitation numbers. Again, e.g., the plateau-structure of 
the doubly excited states confirms the above reasoning. 
These results indicate a second fundamental difference of 
the crystalline structures found in our setup from previ- 
ously proposals. The GSC are located at fixed positions 
relative to the trap borders. In our case, however, the 
resonant excitation crystals are not fixed relative to the 
ensemble, but can float over a certain position range in 
the atom gas from realization to realization. This leads 
to the characteristic differences in the excitation density 
shown in Fig. [31 

To reveal the crystalline character in our setup, a mea- 
surement as illustrated in Fig. 21 could be performed. In 
each run, a certain excitation pattern is obtained. To 
compensate for the fioating, the positions of the leftmost 
excited atoms in the different runs are matched by a shift 
of the position axis. Then, the shifted data of the differ- 
ent runs are averaged. A histogram of 2000 runs analyzed 
in this way is shown in Fig. 01 and clearly shows the crys- 
talline structure at long times. At short excitation times 
the correlation peaks build up successively. Remarkably, 
after 500 ns the qualitative features of the steady-state 
situation are present already. Thereafter, only the abso- 
lute height of the peaks increases. 

It can be seen from Fig.|4]that the correlations decrease 
with increasing distance, unlike in an ideal crystal. But 
the correlations extend over a multiple of the interaction 
length scale, and we have verified for different trap sizes 
up to 50 /im that the crystalline ordering extends over 
the entire ensemble volume for the parameters considered 



in Fig. [TJ In consistence with our interpretations, we 
found that the lattice constant of the floating crystals is 
independent of the ensemble volume. 

Experimentally, laser-excited Rydberg atoms in a 
quasi one-dimensional dipole trap could be used to probe 
our theoretical predictions. As the simplest observable, 
the total number of excited atoms could be measured as 
a function of the detuning and the density. This way, 
the shift of the maximum of the Rydberg excitation to 
higher detunings with increasing density could be verified 
as a first manifestation of the resonant excitations. A 
position-resolved measurement of the excitations would 
allow to verify the predicted resonance peaks in g^'^' , and 
the formation of the crystalline structure. By determin- 
ing the position of the resonances as a function of the 
detuning, also the precise structure of the interaction 
potential could be probed. Another way to access the 
crystalline structure would be a photon correlation mea- 
surement [l7|. Furthermore, the spatial dependence of 
the various excitation-number subspaces arising from the 
fioating nature of our crystals could be analyzed. Finally, 
the resonant excitation crystals are formed dynamically 
throughout the time evolution and do not require a ther- 
malization of the system. This invites a time-resolved 
study of the dynamic formation of spatial correlations. 
Challenges may arise if the trap boundaries are not well- 
defined, leading to a blurring of the predicted position- 
resolved excitation structures. This could be overcome 
by only exciting a well-defined part with near-constant 
density out of the center of a trapped cloud. In approx- 
imately ID systems or higher-dimensional systems, the 
higher-order resonance conditions do not uniquely deter- 
mine the positions of all involved particles. This leads to 
a broadening of the resonances in g^'^\ 

In summary, we have shown that Rydberg excitation 
crystals form at off-resonant laser driving. These crys- 
tals extend over the entire trap volume, and differ from 
previously studied crystals as they are not ground state 
crystals, as their lattice constant is independent of the 
ensemble length, as they are not localized relative to the 
gas, and as the Mandel Q parameter, the Rydberg excita- 
tion density and the total number of Rydberg excitations 
yield qualitatively different results. 
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